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Abstract. The ratio of total mass m* to surface radius r* of spherical perfect fluid 
ball has an upper bound, Gm*/(c 2 r*) < B. Buchdahl obtained Bb uch = 4/9 under 
the assumptions; non-increasing mass density in outward direction, and barotropic 
equation of state. Barraco and Hamity decreased the Buchdahl’s bound to a lower 
value f?BaHa = 3/8 (< 4/9) by adding the dominant energy condition to Buchdahl’s 
assumptions. In this paper, we further decrease the Barraco-Hamity’s bound to 
B new ~ 0.3636403 (< 3/8) by adding the subluminal (slower-than-light) condition 
of sound speed. In our analysis, we solve numerically Tolman-Oppenheimer-Volkoff 
equations, and the mass-to-radius ratio is maximized by variation of mass, radius and 
pressure inside the fluid ball as functions of mass density. 
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A question of self-gravitating compact object is the subject of this paper: How much 
weight or how small radius can a compact object possess? A quantitative answer is given 
by the mass-to-radius ratio Gm */(c 2 r*), where m* and r* are respectively the total mass 
and surface radius of compact object. Assuming static spherical symmetry for simplicity, 
which implies the outside region of compact object is Schwarzschild geometry, this ratio 
is bounded above, Gm*/(c 2 r t ) < B. The upper bound needs to satisfy, B < 1/2 
(<t^ 2 Gm^/c 2 < r*), in order to avoid gravitational collapse. 

An interesting issue of this bound is whether B is less than 1/3 or not. If an 
inequality, B > 1/3 (•<=> 3 Gm^/c 2 > r*), holds for static spherical case, unstable circular 
orbits of photons can appear in the outside Schwarzschild geometry. If such super¬ 
compact object, which possesses unstable circular orbits of photons but no black hole 
horizon, do neither emit nor reflect any radiation, then we cannot distinguish it from 
black holes by observing the so-called black hole shadows. Here, the black hole shadow 
is a dark region which is expected to appear in a fine image of optical/radio observation 
of black holes (see BE] and references therein): Consider a case that an optical source 
is extended behind a black hole and does not enter the inside of unstable circular orbits 
of photons. Observer cannot detect photons which passed through unstable circular 
orbits inward, because those photons are absorbed by black hole eventually. Hence, 
in the image of optical source, there should appear the dark region, on which those 
photons would be detected if black hole did not exist. This dark region is the black hole 
shadow, and the boundary of shadow is determined by the photons propagating on null 
geodesics winding many times around unstable circular orbits of photons. This means 
that a super-compact object possessing unstable circular orbits but no black hole horizon 
can provide us the same optical image of shadow with black holes. Since the resolution 
of image by radio observation is now approaching the visible angular size of largest 
black hole candidate M, an investigation of possible super-compact object seems to 
be an important issue for near future observational study of black holes. Therefore, an 
interesting and important issue is whether B < 1/3 holds or not. 

Some exotic models of super-compact objects have been proposed such as 
gravastars, boson stars and so on. The gravastar as super-compact object has already 
been examined [5j, whereas the others remain to be examined. Those exotic models 
may be interesting. However, in this paper, we focus on a rather usual model. 

Assuming that (i) the compact object is a static spherical ball of perfect fluid, (ii) its 
mass density is non-increasing with respect to radial coordinate, and (iii) its equation of 
state is barotropic, Buchdahl obtained Gm*/(c 2 r t ) < Bb uc h = 4/9 [6, 7j. Furthermore, 
by adding the dominant energy condition to Buchdahl’s assumptions, Barraco and 
Hamity decreased the Buchdahl’s bound to a lower value £?B a Ha = 3/8 (< 4/9) pj. 
However, this upper bound is greater than 1/3 (< £>B a Ha)- 

Other works on the bound B have been performed. For example, an effect of 
cosmological constant is examined [9], a case with infinite surface radius is analyzed [10], 
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and a tangential pressure is considered m However, the Barraco-Hamity’s bound, 
^BaHa = 3/8, has not been lowered so far. 

As explained in detail in section [2l we add the subluminal (slower-than-light) 
condition of sound speed to the assumptions of Buchdahl and Barraco-Hamity. We 
discuss, in section [21 how the subluminal-sound-speed condition restricts the form of 
equation of state, and also that this condition is more restrictive than the dominant 
energy condition. Then, our upper bound is lower than the Barraco-Hamity’s bound, 

Cirirn / S \ 

-T 1 < I3 new ~ 0.3636403 (< B baHa = x) • (1) 

The reason why £> new is approximate value is that we have performed numerical 
integration of Tolman-Oppenhcimer-Volkoff equations. Although the upper bound is 
lowered, it still remains greater than 1/3 (< £> new ). Therefore, at present, we cannot 
deny the possibility that the same shadow image can be obtained from black holes and 
some super-compact object made of fluid matter. In order to sharpen the upper bound 
of Gm»/(c 2 r»), some other condition should be assigned to the form of equation of state. 
(We will report the case of polytropic equation of state in the other paper.) 

Section [2] is devoted to the details of our analysis. Section [3] is for conclusion and 
discussions. 


2. Variational analysis of mass-to-radius ratio 


2.1. Definitions and assumptions 


As explained in section [T] we are interested in a static and spherically symmetric perfect 
fluid ball. A line element of spacetime is 

dr 2 

ds 2 = g tw dx> 1 dx" = -e^Vdf 2 + 1 — . ... _ . + r 2 (d9 2 + sin 2 Odip 2 ) , (2) 

1 — 2Gm(r)/(dr) 

where (t,r,9,(p) is spherical poler coordinates, <F(r) gives a lapse function, and m(r) is 
a mass of perfect fluid contained in spherical region of radius r. The stress-energy- 
momentum tensor of perfect fluid is T^ u = a(r)c 2 u,j,u u + p{r) {g, JV + u^Uv), where 
u = e~^d ct is a four-velocity of static perfect fluid, and cr(r) and p{r) are respectively a 
mass density and pressure of perfect fluid. 

By the regularity of spacetime at centre, a condition m( 0) = 0 holds. This implies 
that the mass density at centre cr c = <r(0) is finite. We normalize all quantities by a ci 


„ vux 

R := -r 


S(i?) := 


a(r) 


M(R ) := 


V& 


cr r 


m(r) 


P(R ) == 


p(r) 


( 3 ) 


<7 C C° <J c C* 

These are dimension-less. The lapse function, <h(r) := $(/?), does not need 

normalization because <f> is originally dimension-less by definition (J2]l. 

Following Buchdahl 01Z], we adopt two assumptions. One of them is a barotropic 
equation of state, 


P = P(E). 


( 4 ) 
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Hereafter, we regard the mass density £ as an independent variable, and the others are 
functions of it, 


R = R{E) , M = M (£) , P — P(£) , $ = $(£). (5) 


Another assumption is a non-increasing mass density in outward direction, 
dP(£) 


d£ 


< 0 . 


( 6 ) 


This implies the correspondence between variables R and £ is one to one. 

The surface of fluid ball is defined by vanishing pressure, where the hydrostatic 
equilibrium holds between fluid ball and outside vacuum region. Then, the mass density 
at surface £* is determined by 


P(£*) — 0, 


(7) 


and the total mass M* and surface radius P* of fluid ball are respectively given by 


M* = M(£*) , P* = P(£*). (8) 

At the centre of fluid ball, the mass density is unity £ c = 1, and conditions, P(l) = 0 
and M( 1) = 0, should hold. Note that the mass density takes values in an interval, 

£* < £ < 1, (9) 


where the surface mass density satisfies 0 < £* < 1. 

The outside region of fluid ball, P > P*, is Schwarzschild geometry of mass M*. 
The inside region, P < P*, is determined by the Einstein equation and conservation law 
T ,w . u = 0, which are reduced to Tolman-Oppenhcimer-Volkoff (TOY) equations, 


d M(£) 
d£ 

d P(S) 

d£ 

d$(£) 

d£ 


4vrP(£) 2 £ 


dP(£) 

d£ 


A{M, P, P; £) 
1 


dP(£) 


d£ 
d P(£) 


£ + P(£) d£ 


(10 a) 
(106) 
(10c) 


where 


A{M,R,P]E) : = 


; £ + P(£) ] [ M(£) + 4vrP(£) 3 P(£) ] 


( 11 ) 


P(£) [P(£) — 2M(£) ] 

Given a concrete functional form of equation of state, two functions P(£) and M(£) 
are obtained by solving (llOd) and (11061) . Substituting those solutions into (llQcl) . $(£) 
is obtained. 

In addition to Buchdahl’s assumptions (J4J) and (J6]), we assume subluminal condition 
of sound speed, 

dP(£) 


d£ 


< 1 . 


( 12 ) 


Given the above formulation, the remaining freedom is the functional form of P(£) 
under the condition (1T2|) . Since the solution of TOV equations depend on the concrete 
form of P(£), the value of total mass M* and surface radius P* of fluid ball vary with the 
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concrete form of P(E). In following sections, by the variation of functional form of P(E) 
under the condition (1121) . a possible upper bound B of the ratio M*/P* (= Gm*/(c 2 r*)) 
will be calculated, 


M* „ 

-rT<B. 


(13) 


Using such variational method, our aim is to obtain the value of B lower than the 
Barraco-Hamity’s bound PeaHa = 3/8 which is lower than the famous Buchdahl’s bound 
^Buch = 4/9. 

Here note that Barraco and Hamity [8] considered not only the Buchdahl’s 
assumptions (HJ) and (Jf)J) but also the dominant energy condition which gives the following 
inequality for the perfect fluid ball, 


P(E) <E. (14) 

We expect that some physically reasonable (non-exotic) condition, which is more 
restrictive than the dominant energy condition, can decrease further the Barraco- 
Hamity’s bound to a lower value. Concerning this expectation, let us point out that 
the subluminal-sound-speed condition and definition of surface (J7J) predict the 
inequality (TT4|) as well. The subluminal-sound-speed condition restricts not only the 
value of pressure such as inequality (TT4l) but also the value of its differential such 
as inequality (fTIj) . This implies that the subluminal-sound-speed condition is 
more restrictive than the dominant energy condition. Therefore, we expect that the 
subluminal-sound-speed condition is sufficient for us in search of an upper bound lower 
than Barraco-Hamity’s bound. 


2.2. Variational method: Strategy 


The variational method which we are going to use is known as the optimal control 
theory. It has already been applied to an estimation of maximum mass of neutron star 
by Rhoades and Ruffini p21J. (See also (13[ E] for more details about application of 
the optimal control theory to neutron star mass.) We apply the optimal control theory 
to the search of upper bound B of inequality (fT3l) . However, we do not require readers 
to have knowledge of optimal control theory. All ideas of the variational analysis are 
explained below. 

We introduce an auxiliary variable [/(E) by 

sin 2 U(S) := , ( 15 ) 

Then, the subluminal-sound-speed condition (TT2l) is automatically satisfied, sin 2 U < 1. 
In order to make use of variational method, it is useful to express M*/P* by an integral 
form, [jj 


M* 

R * 



d (M/R) 
d R 



( d {M/R) \ 
V dE ) 



(16a) 


f The signature of (I16al) is not — f dRd{M / R) / dR but + / dRd(M/R)/dR, because of an inequality 
d(M/R)/dR > 0 which is found by M/R ~ R 2 . 
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where the integrand L is arranged to the following form by using relation (1131) and TOV 
equations (1 10 a|) and (1106|) . 

LMR'P.UM- (164) 

Our problem is to maximize the functional (1 16 all under three constraints; the relation 
(TT5j) and TOV equations (IIPop and (IK) 61) . In following discussion, (1 10 c|) is not considered, 
because $ does not appear in the integrand L. 

Because there are three constraints, we use the Lagrange’s multiplier-method. 
Define a functional, 

/:= /dS[L + y M (S) Cm + V R (E) Or + y P (E) C P ], (17a) 

as* 

where functions 1 m(S), V r (E) and 1 P (E) are Lagrange multiplier, and Cm, C r and C P 
are defined according to the constraints fllPaj) . (11061) and (fT5j) . 


C m (M,R,P; E) = 
C r (M,P,P;E) = 


4ttP(E) 2 E) . 2 dM(E) 


A(M, R, P; E) 


A(M, R, P; E) 


sin 2 67 (E) - 
sin 2 67 (E) 


dE 

dP(E) 

dE 


C P (M, R, P; E) = sin 2 f/(E) 


dP(E) 

dE 


(176) 
(17c) 
(17 d) 


Extremal value of is given by the condition 51 = 0 under variations, 

M(E) M(E) + 5M(E) 

R{ E) ->■ M(E) + 5i2(E) 

P(E) -> M(E) + 5P(Yi) 

67(E) ->• 67(E) + 567(E). 

Given the functional expression I of M*//?,*, the strategy of our variational analysis 
consists of following steps: 


(18) 


stepl Divide the interval (J9]) of E into many infinitesimal intervals. Then, solve the 
equation 51 = 0 in order to maximize the functional I at each infinitesimal interval. 
step2 Integrate the maximized / of every infinitesimal interval of E. Then, we 
will obtain a formal expression of the upper bound B > M*/P* for the global 
interval (|9l) . 

step3 Calculate numerically the formal expression of B obtained in step2. 


2.3. Variational method: Stepl 

We divide the interval (JHD into many infinitesimal intervals. Look at one infinitesimal 
interval of E, 

S d <E<E u , (19) 

where E* < Ed < E u < 1 and E u — Ed -C 1. Then, in the stepl of our variational 
analysis, we maximize by solving 51 = 0 under the variations (fT5D in the 
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Figure 1. Extremization of in an infinitesimal interval Ed < E < E u . An 

original equation of state expressed by dashed curve is deformed to an extremal case 
of which forms an parallelogram. 


interval f[T9|) . Note that, since we are considering a single interval (1T9|l in the stepl, 
we fix variables (M, R , P, U) in the remaining intervals, E* < E < E d and E u < E < 1. 
This indicates boundary conditions, 

4X(E d ) = 0 , <5X(E U ) = 0 (X = M, R, P, U) . (20) 


The Enler-Lagrange equations of 51 = 0 are followings: Variation of U gives 
-M 


-47Tfl s E ,. 4jrfl 2 E ,, 1 ,. 

+ u ~— + R I + P 


sin U cos U = 0 , 


( 21 ) 


variations of X(= M, R , P) give 

<9>x(£) 


as 


d[(M-4nR 3 Y)/(AR 2 )] d[AnR 2 Y/A) 0(1 /A) 

-1- AM--b Yr- 


sin 2 U , (22) 


dx ox ' ox 

and variations of Yx give Cx = 0 which are constraints; TOV equations (llOap . (110 61) . 
and relation (TT3T) . Equation (0T]1 gives 


U = 0, -, 
’ 2 ’ 


(23) 


or an algebraic equation given by vanishing the inside of square bracket of m • 
However, the latter case (vanishing square bracket) is impossible under the boundary 
condition (12Uj) as shown in Appendix A Therefore, U should be constant given in 
which denotes for extremal case of 


d PP) 

dE 


0 , 1 . 


(24) 


This gives a constant pressure or linear equation of state, P(E) = E — E*. It is already 
revealed by Nilsson and Uggla |15] that, for the linear equation of state (P = E — E*), 
the total mass M* and surface radius P* are both finite for E* ^ 0, but both infinite for 
E* = 0. However, the behaviour of M*/P* has not been investigated so far. 


§ In the case of Buchdahl’s analysis [B], the maximum value Af* //?,* = t?Buch = 4/9 is realized for the 
fluid ball of constant mass density. The constant mass density implies that the sound speed is infinity. 
On the other hand, we are considering the subluminal-sound-speed condition. Therefore, the value 1 
of sound speed for the extremal case of M^/R^ in our analysis is consistent with Buchdahl’s analysis. 
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Due to the extremal condition (1241) and boundary condition (12U|) . a parallelogram 
appears in S-P plane as show in hgure [0 The equation of state given by edge 
(a) + (b) or (c) + (d) corresponds to a maximum or minimum value of M*/P* under 
the variations (TT8|) in the infinitesimal interval ([T9]) . In order to judge which edge 
corresponds to maximum or minimum, we calculate their difference, 

:= [ dS L- f dEL, (25a) 

V R * ' J (c)+(d) J( a)+(b) 

where the hrst term corresponds to an integral with equation of state given by the edge 
(c) + (d), and the second term is an integral with the edge (a) + (b). Here, note that 
TOV equations (110 al) and (I106j) result in that M (E) and P(E) are constant for U = 0 
(dP/dE = 0). Therefore, the edges (b) and (c) cause no contribution to M*/P*. Then, 
the difference A (M*/P*) is calculated to be 


A 





L(M, R,P,U = 1; E d + AS) 5 S - L(M, R, P, U 
dL(M , R,P,U = l;E d ) 


<9E d 

P d — 2 M d 
P d (E d + P d ) 2 


AS AS 


hS AS > 0 , 


l;S d )5S 


(256) 


where 6S and AS correspond to two edges of parallelogram as shown in hgure [T] and 
M d , P d and P d are mass, radius and pressure at E d . This inequality (12561) denotes that 
the equation of state with edge (c) + (d) (edge (a) + (b)) corresponds to the maximum 
(minimum) of M*/P* under variations in the infinitesimal interval (fT9l) and boundary 
condition (12U|) . 


2. f. Variational method: Step2 


Next, we proceed to an analysis in the global interval ([2]) of S. Consider a case that 
the values of pressure at centre, P c , and mass density at surface, E*, are fixed. In this 
case, we can repeat the deformation of equation of state using an infinitesimally small 
parallelogram. Finally, the functional form of equation of state in interval ([9]) reaches 
one of following two options: 


P( s) = 


S — E* in E* < S < E* + P c 
P c in E* + P c < S < 1 


P(£) = 


0 in E* < E < 1 — P c 

E — (1 — P c ) in 1 — P c < E < 1 


(26 a) 
(266) 


These forms are shown in hgure [21 No contribution to M*/P* causes from parts, P = 
constant, as explained before calculating (125 61) . The inequality (125 6j) for infinitesimally 
small parallelogram indicates that the equation of state (12661) gives the maximum value 
of M*/P* for given P c and E*, whereas (126 al) gives the minimum for given P c and E*. 
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Figure 2. Maximization and minimization of A/*/it* in the full interval E* < E < 1 
with fixing P c and E*. Horizontal dashed parts of e.o.s. do not contribute to Af*/i?*. 


We should emphasize that the above discussion is applicable for only the case with 
fixed P c and £*. It has not been examined so far whether or not M*/P* for equation 
of state (1266)) (or (126a)) ) is maximum (or minimum) even when the values of P c and £* 
vary. There is a possibility that a value of M*/P* for equations of state (1 2 Gal) or (12661) 
for certain values of P c and £* is neither maximum nor minimum under variation of P c 
and £*. 


In order to find a true maximum value of M*/P* under variation of P c and £*, we 
regard M*/P* for equation of state (1266)) as a function of P c , 


/(-Pc) : = 


M* 
R * 


for equation of state (1266)) . 


(27) 


Here, note that the equation of state (12661) in the interval, 1 — P c < £ < 1, depends 
on only P c . The maximum value of f(P c ), which is denoted by B, is the desired upper 
bound of mass-to-radius ratio, M*/P* < B. A numerical plot of /(P c ) is going to be 
shown in next section fUA 

In order to find a true minimum value of M*/P* under variation of P c and £*, we 
regard M*/P* for equation of state (126 ah as a function of P c and £*. The minimum 
value of this function is the lower bound of M*/P*. Here, note that this lower bound can 
be read from f(P c ) by following discussion: For the equation of state ([26 al) for given P c 
and £*, the quantity £* + P c can be regarded as a central mass density!]]) Then, following 
the normalization (]3|), we transform variables by £ = £/(£* + P c ), P = P/(£* + P c ), 
M = \/£* + P C M and R = \/£* + P c R. This transformation changes the form of 
equation of state (126 qj) in interval, £*<£<£* + P c , to the following form, 


P(£) = £ - (1 - P c ) for 1 - P c < £ < 1 


(28) 


where P c := P c /(£* + P c ). Furthermore, the form of TOV equations for transformed 
variables are the same with QlOaj) and (1106)) . Hence, the value of M*/P* for equation of 
state (l28l) is equal to the value of M*/P* for equation of state (126 qj) . On the other hand, 
the form of equation of state (128|) is the same with equation of state (126 6[) in interval, 


| We should also remember, R( E) = constant for P = constant due to TOV equations. 
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Figure 3. M*/P* as a function of P c for linear equation of state (12661) . Differential is 
plotted by an approximation, d/(P c )/dP c ~ [f(P c + S) — f(P c —S)]/S, where <5 = 10 -6 . 
These plots are made with Mathematica ver.10. 


1 — P c < E < 1. Therefore, the value of M*/ R t for equation of state (126 all for given P c 
and £* is equal to that for equation of state (126611 with replacing P c by P c . This implies 
that the value of for equation of state (126 aft can be read from f(P c ), and the 

minimum value of f(P c ) is the lower bound of 


2.5. Variational method: StepS 


The function f(P c ), defined in (l27jl . can be calculated by following steps: 


(i) Solve numerically TOY equations (1 10 all and (110611 with the equation of state (126611 


for given value of P c . A technical remark is summarized in Appendix B 


(ii) Calculate from numerical solutions of step (i). 

(iii) The function f(P c ) is calculated by repeating steps (i) and (ii) for different values 
of P c in the interval, 0 < P c < 1. 


Our numerical result is shown in figure [3] A smooth graph of f(P c ) is obtained. From 
the plot of differential of f(P c ) in figure [3j we find that the maximum of f(P c ) is given 
at 


Pc = Abound ^ 0.8386058 . (29) 

Then, as discussed in previous section l2~4l the upper bound B new > M^/R* is given 
by the maximum value, B ncw ~ /(-Pbound), which is evaluated numerically to be 
£> new — 0.3636403. This upper bound gives our conclusion (PQ). On the other hand, 
the lower bound read from figure [3] is zero given at P c = 0. No finite lower bound is 
found in our analysis. 
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By adding the subluminal-sound-speed condition (fT2l) to Buchdahl’s assumptions, we 
obtained figure [3] and the bound of mass-to-radius ratio of a barotropic fluid ball, 


0 < 


-M* 

P* 


< B nPW ~ 0.3636403 . 


(30) 


This upper bound is lower than Barraco-Hamity’s bound £> Ba Ha = 3/8. However, our 
upper bound is greater than 1/3 (< B new ), and there remains a possibility of an existence 
of super-compact object possessing unstable circular orbits of photons but no black hole 
horizon. 

As read from figure [3j the upper bound B new corresponds to the central pressure, 
P c = Pbound °f O- This value of central pressure may be understood as a result of 
trade-off between pressure’s two effects; a contribution to mass density (attractive force) 
and an effect pushing fluid outward (repulsive force). 

From the above, the maximum value of M*/P* is B new of (1301) which is obtained 
by solving the TOV equations with the central pressure P c = Pbound of (1291) and linear 
equation of state (12661) . Furthermore, by our numerical integration of TOV equations, 
the mass and radius of fluid ball in the case of maximum mass-to-radius ratio are 

Af* (Pbound) = Afbound — 0.2014578 
R* (Pbound) = Pbound — 0.5540028, 

which gives our bound £> new = Mb OU nd/Pbound- From these values together with the 
normalization (J3J), we can estimate typical un-normalized values of physical quantitis 
for the case of maximum mass-to-radius ratio. With denoting the un-normalized mass 
in unit of solar mass, rebound = ckMq, where a is a dimension-less factor, then the 
un-normalized central mass density abound, central pressure Abound and radius r boU nd are 
calculated as 


(31) 


Abound 


C 6M bound 

Ground 


Pbound Abound C -Abound 

_ G^bound 

Abound 


~ x 10 16 g/cm 3 

cr 

~ x 10 39 GeV/cm 

cr 

~ 4.1 a km. 


(32) 


n -2 

' ^new 

Note that, the case of maximum mass-to-radius ratio corresponds to the linear 
equation of state (126 6j) giving the luminal (equal-to-light) sound speed inside the fluid 
ball. If the equation of state, P(£), will be restricted to some functional form other 
than the linear form, then the upper bound B of M*/P* will take a value different from 
£>new of (1301) . We will report the case of polytropic equation of state in the other paper 0 

We will show in the other paper that, when the equation of state is restricted to be the polytropic one 
(P(£) oc £ 1+1 / ra ) under the sublunrinal-sound-speed condition, the maximum value of M*/R* becomes 
a value lower than B new of (H(l . M^/R* < 0.281. In deriving this result of polytrope, we will not use 
the variational method, but use the other numerical search method. 
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Finally we make a mathematical comment on a limiting behaviour found from 
figure El 

lim f(P c ) = j . (33) 

P c ->1 4 

At the limit P c —> 1, the equation of state (126 M asymptotes to a simple form P(£) = £ 
(0 < £ < 1). Further, functional forms of M(£) and P(£) converge to an exact solution 
of TOY equations with equation of state P(£) = £ ; 


M(E) = ifl(E) = . 


(34) 


This results in a finite limit (133]h However, the total mass and surface radius of fluid 
ball diverge M* —> oo and P* —* oo as P c —> 1, because the surface mass density is zero, 
£* —» 0 as P c —> 1. This result is consistent with Nilsson-Uggla’s numerical result [15] 
that, for the linear equation of state, P = £ — £* (£* = 1 —P c ), the mass M* and radius 
P* are both finite for £* ^ 0 but diverge as £* —> 0. 
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Appendix A. Impossibility of vanishing square bracket of (12111 


In this appendix, we show that, as a solution of (]2Tlh the vanishing square bracket of 
(l2Tf is impossible under the boundary condition (120]) . 

Assume that a solution of ()2TTl is given by vanishing the square bracket, we find 

1 '• - -- ‘ • (A.l) 


rp A(M, R, P; E) 

On the other hand, we obtain from 

dlM R + 8tt R 3 P 


r m 9 

/ 1 \ 1 

Y R + — + 4nR 2 X 

1 

rOh 

1_ 


d£ 

din 

d£ 


dip 


-Y P sin 2 U 


sin 2 U 


(P- 2M)(M + 4vrP 3 P) 

1 1 12ttP 2 P \ 

“ ' P + (P^2M) _ M + RHT ) 


2 M 


4tt£ (l - 2PY 


M 


A(M, P, P; £)P 2 
Y P sin 2 U 
sin 2 U 


P 3 

M£ + 4ttP 3 (P£ + £ + P) 


A(M, P, P; £) 
lp sin 2 U. 


(A.2) 


d£ (£ + P)(M + 4 ttP 3 P) 

Equation (IA.4|1 with boundary condition (l20j) results in lp = 0. Then (jA.ljl gives 


(A.3) 

(A.4) 
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Differentiate this equation by £ and substitute (I A. 21) . (I A. 31) and TOV equations (110 al) 
and (110 61) . we obtain 

1 ,, M(£) 


>m(£) = 


Yr&) = 


(A. 6 ) 


P(£) ’ ' P (£) 2 ’ 

where 1 r is obtained by substituting Ym in (IA.5H . These solutions cannot satisfy the 
boundary condition ([20]) . Hence, the solution of (l2Tj) is given by 


Appendix B. On numerical treatment of TOV equations 


Right-hand sides of TOV equations (110 q j) and (I106|) are indeterminate form at centre 
because of the conditions, M —> 0 and R —> 0 as £ —» 1. Therefore, in solving 
TOV equations numerically in step (i) of section [2751 we have made use of perturbative 
solutions near the centre. 

In order to consider a perturbation near centre, we regard the radius R as an 
independent variable, and the mass density as a function of radius, £(P). TOV 
equations (llOal) and (I106j) are rearranged to the following forms, 


d M(R) 


= 4vr P 2 £(P) 


dP(R) [ £(P) + P(R)} [ M(R) + 4t tR 3 P(R) ] 


(B.l) 


(B.2) 


d R ’ d R R[R-2M(R)} 

For a sufficiently small radius R 1, we introduce perturbations, 

M(R ) = M (1) R + M (2) i? 2 + M (3) P 3 + • • • 

P(R ) = P c + P(i)R + P{2 )R 2 + P(3 )R 3 + • • • 

£(i?) = 1 + £(i)i? + £(2 )R~ + £(3 )i ? 3 + • • • , 

where conditions M(R = 0) = 0, £(i? = 0) = 1 and P(R = 0) = P c are included. 
Substituting (IB.2D into (IB.ID . we obtain Mp) = 0, M 2 = 0, Pp) = 0 and remaining 
parts, 

M{R) = ~7rP 3 + 7t£ (1 ) P 4 + • • • 
o 


P(R) - P, - + 3P c ) (1 + P c )Jt 2 - y<7 + 15P C )E (1 ,B 3 + • • ■ < B - 3 ) 


A(1 + 3 P 0 ) (1 + f 

£(P) = 1 + £p)P + £( 2 )P^ + £(3 )P 3 + • • • , 

where the central pressure P c and coefficients £( n ) (n = 1, 2, 3, • • •) are determined by 
the concrete form of equation of state. 

Comparing P(R) in (IB. 3D with equations of state (12661) . we obtain £p ; = 0 and 
£( 2 ) = P ( 2 ) • Hence, denoting a small radius by R$ -C 1, the mass density £5 and mass 
Ms at R = Rs are approximately given by 


£,5 = 1 - -7r(l + 3P C ) (1 + P C )R] , M d = -ttR 3 s 


If the mass density near centre £5 is given, the others are determined by 

3/2 

/ .til - / .I I *-+ / .til - / jJT I \ 

Rs = 


3(1-£* 


27t( 1 + 3P C )(1 + P c ) 


3(1 — £ 5 ) 


4 

Ms = — 7T . 

3 V2tt(1 + 3P c )(1 + P c ) 


(B.4) 


(B.5) 
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In step (i) of section 12.51 we have solved TOV equations (11 Gap and (11061) in an interval, 
£* < £ < T,$, with initial condition (IB.5|) . Also, the convergence of numerical solutions 
has been checked by varying £ 5 . Our results in figure [3] are made using £5 = 1 — 10 -6 . 
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